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Abstract: NOz^A experiment has reoptimized its event selection criteria in light of the 
recently measured moderately large value of #13 . We study the improvement in the sensitiv- 
ity to the neutrino mass hierarchy and to leptonic CP violation due to these new features. 
For favourable values of Sep, NOz^A sensitivity to mass hierarchy and leptonic CP violation 
is increased by 20%. Addition of 5 years of neutrino data from T2K to NOz^A more than 
doubles the range of Sep for which the leptonic CP violation can be discovered, compared 
to stand alone NOz^A. But for unfavourable values of 5qp, the combination of NOz^A and 
T2K are not enough to provide even a 90% C.L. hint of hierarchy discovery. Therefore, 
we further explore the improvement in the hierarchy and CP violation sensitivities due to 
the addition of a 10 kt liquid argon detector placed close to NOz^A site. The capabilities 
of such a detector are equivalent to those of NOz^A in all respects. We find that combined 
data from 10 kt liquid argon detector (3 years of v + 3 years of v run), NOz^A (6 years of 
v + 6 years of v run) and T2K (5 years of v run) can give a close to 2a hint of hierarchy 
discovery for all values of <5cp- With this combined data, we can achieve CP violation 
discovery at 95% C.L. for roughly 60% values of 5qp- 
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1 Introduction 

The recent measurements of #13 by the reactor neutrino experiments Double Chooz [1, 2], 
Daya Bay [3, 4], and RENO [5, 6] are a very welcome news for future neutrino oscillation 
experiments. In addition to confirming the non-zero #13 as hinted earlier by the accelerator 
experiments T2K [7, 8] and MINOS [9, 10], these experiments have determined the value 
of #13 quite precisely. Global fits to the world neutrino data [11, 12] give a best fit value 
sin 2 2#i3 = 0.095 with a la uncertainty of about 10% [12]. Daya Bay experiment is expected 
to reduce the uncertainty to 5% level by the time it finishes running in 2016 [13]. Since 
the value of sin 2 26is is moderately large (in fact, it is just below the CHOOZ upper 
limit [14, 15]), the current (T2K) [16] and upcoming (NOz^A) [17] experiments now have a 
chance of determining the remaining unknowns of neutrino oscillations: i.e. (a) neutrino 
mass hierarchy (equivalently the sign of Am^ = m 2 — m 2 ) and (b) existence of CP violation 
in the leptonic sector. 

CP violation in the leptonic sector has drawn tremendous interest because of the possi- 
bility of leptogenesis leading to baryogenesis and the baryon asymmetry of the universe [18]. 
Leptogenesis requires the existence of CP violation in the leptonic sector; see [19] for a re- 
cent review. The possible connection between leptogenesis and neutrino oscillations has 
been discussed in [20-22]. It is likely that the CP violating phase in neutrino oscillations 
is not directly related to the CP violation leading to leptogenesis. But a demonstration 
of CP violation in neutrino oscillations provides a crucial guidepost for models of leptonic 
CP violation and leptogenesis. 
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The experiments NOz^A and T2K expect to achieve the above goals by measuring the 
- > v e oscillation probability and its charge conjugate — > v e oscillation probability 
Pfie- Since these experiments have moderately long baselines (295 km for T2K and 810 km 
for NOz/A), the matter effects due to neutrino propagation through Earth [23] are impor- 
tant, especially for NOz/A [24]. The matter term modifies P^ e (and also Pp,g) differently 
for normal hierarchy (NH where Am^ is positive) and for inverted hierarchy (IH where 
Anij] is negative). Thus these experiments are capable of distinguishing between the two 
hierarchies. The matter term also changes sign when we switch from neutrino mode to 
anti-neutrino mode. Hence the matter effects induce a CP like change in the oscillation 
probabilities. Thus we have an entanglement of changes caused by matter effects and the 
genuine CP violating phase 5cp- Due to this, we get two degenerate solutions for a given 
experiment: one with the right hierarchy and the right value of <5cp and one with the wrong 
hierarchy and a wrong value of 5cp- These degenerate solutions can be unravelled if data 
from two different experiments with different baselines is available [25-30]. 

It was shown that the data from T2K and NOz/A should be synergistically combined 
to obtain the best possible sensitivity to hierarchy [31]. The presently planned runs of T2K 
and NOz/A will not be able to determine the hierarchy for the whole range of <5cp [32-34]. 
Matter effects increase P(va —> v e) for NH and decrease it for IH and vice verse for Pn g . 
For 5qp in the lower half-plane (LHP, —180° < <5cp < 0), Piy^ — > v e ) is larger and for 
<5cp in the upper half-plane (UHP, < 5qp < 180°), P{v^ — > u e ) is smaller. Hence, for 
the combination (NH, LHP), the values of P(vn -4 z/ e ) are much higher than those for IH 
(and Ppg values are much lower). Similarly, for the combination (IH, UHP), the values of 
P{ v u ~4 v e) are much lower than those of NH (and Ppg values are much higher). Thus, 
LHP is the favourable half-plane for NH and UHP is for IH [35]. Therefore, NOz/A by itself 
can determine the hierarchy if 5cp happens to be in the favourable half-plane [33, 35]. The 
range of <5cp values, for which the hierarchy can be determined, increases to some extent 
if the T2K data is included. It was shown that the hierarchy can be determined for the 
entire range of 5cp if there is 50% more data from NOz/A and twice the data from T2K 
[35]. 

In light of the moderately large value of #13, NOz/A experiment has reoptimized their 
event selection criteria, with more events in both signal and background [36, 37]. These 
new criteria improve the hierarchy determination ability significantly in the favourable 
half-planes of <5cp but not in the unfavourable half-planes. The hierarchy sensitivity, for 
the unfavourable <5cp half-planes, can be improved only with a much larger data sample. In 
this paper, we explore the expected improvement in the sensitivity to both hierarchy and 
leptonic CP violation because of increased data from experiments with 810 km baseline. 
Liquid argon time projection chamber (LArTPC) are emerging as a very good option for 
neutrino detector technology for future neutrino experiments both in Europe [38-41] and 
in the USA [42-44]. Recently, LBNE collaboration has explored the possibility of building 
a 10 to 30 kt LArTPC at either NOz/A or MINOS site [45]. We calculate the combined 
sensitivity of a 10 kt LArTPC at the NOz/A site, in conjunction with NOz/A and T2K 
detectors, to hierarchy and to the existence of leptonic CP violation. 

The paper is organized as follows. In section 2, we discuss the detectors we consider in 
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this analysis, paying particular attention to the improvements implemented in NO^A [36] 
and contrasting them with old NO^A event selection. We also detail how these improve- 
ments have been implemented in our simulations. We briefly describe the LArTPC and 
compare its properties with those of NOz^A. Sections 3 and 4 describe the event spectrum 
and our numerical procedure respectively. In section 5, we present our sensitivity results for 
the determination of hierarchy and the detection of CP violation for various combinations 
of experiments. Finally, in section 6, we mention our conclusions. 

2 Detectors 

In this section, we briefly describe the main features of the detectors considered in this 
report. 

2.1 NOuA 

NOz^A is a 14 kt totally active scintillator detector (TASD) at a distance of 810 km from 
Fermilab, at a location which is 0.8° off-axis from the NuMI beam. Because of the off-axis 
location, the flux of the neutrinos is reduced but is sharply peaked around 2 GeV. This 
leads to two important advantages: 

• The peak flux is close to the first oscillation maximum energy of 1.7 GeV. This leads 
to a large number of signal events. 

• The most problematic background is neutral current (NC) interactions which mostly 
consists of the single ir° production. However, the measured energy of this background 
is shifted to values of energy below the region where the flux is significant. Hence 
this background can be rejected using a simple kinematic cut. 

The experiment is scheduled to run for three years in neutrino mode and three years in 
anti-neutrino mode with a NuMI beam power of 0.7 MW, corresponding to 6 x 10 20 protons 
on target per year. 

Previously, the event selection criteria were optimized to have the least background 
even under the most pessimistic assumption of sin 2 2#i3 = 0. This gave the best possi- 
ble background rejection, at the cost of reduced signal efficiency, which was required in 
the event of very small sin 2 26*13. With the present, moderately large value of sin 2 2#i3, 
NOz^A has relaxed the cuts for the event selection criteria which allow more signal events 
along with more background events. Additional backgrounds, mostly NC, are reconstructed 
at energies lower compared to the true neutrino energy and can be managed by a kine- 
matic cut. The main differences between the old criteria and new criteria are listed below 
[34, 36, 37]. 

• The signal efficiencies for new NOzaA. are higher than that of old one by roughly a 
factor of 2 for neutrino events. We now assume 45% signal efficiency for both neutrino 
and anti-neutrino events as opposed to 26% for v and 40% for v previously. 
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• The background acceptance has also increased. For the NC interactions, present value 
is about 7 times (2% vs. 0.3%) for v and 3 times (3% vs. 0.9%) for u, compared to 
the old criteria. For misidentified muons also, it is about 6 times (0.83% vs. 0.13%) 
and 2 times (0.22% vs. 0.13%) the older numbers, for v and v respectively. 

• Earlier, the number of NC background events was moderate. In such a situation, 
using a Gaussian energy resolution function to obtain the smearing of the background 
events is not a bad approximation. At present though, the NC backgrounds are higher 
by a factor of 5 and but their measured energy, in general, will be in a range below the 
region of large flux. The NC spectrum shift to the measured energies is implemented 
through migration matrices. 

• In our analysis, we have taken care of the neutrino contamination in the anti-neutrino 
beam for both appearance and disappearance channels. It should be stressed that 
while anti-neutrino contamination can be ignored in the neutrino beam, the reverse 
is not true. 

The above optimization criteria were developed by using the event spectra for the 
case of — > v e vacuum oscillations with <5cp = and maximizing the signal events while 
keeping the background events relatively small. We calculated the same event spectrum 
under the same assumptions and adjusted the signal and background efficiencies in the 
simulations, until we obtained the same number of signal events and background events as 
in [36]. 

In Table-1, we summarize the main characteristics of TASD and LArTPC which we 
discuss next. 



Detector Characteristic 


LArTPC @ NO// A 


TASD @ NO^A 


Fiducial Mass 


5 or 10 kt 


14 kt 


Neutrino energy threshold 


0.5 GeV 


0.5 GeV 


Detection efficiency 


85% for m ± (CC), 80% for e± 


100% for /z^CCQE), 45% for e ± 


Energy resolution (GeV) 


0.1^/E/GeV for CC and e ± 


O.Oe^E/GeV for CC ^ and 




sample 


0.085^E/GeV for e ± sample 


NC backround smearing 


Migration Matrices 


Migration Matrices 


NC background acceptance 


1% (both v and v) 


2% (y), 3% (y) 


Mis-ID muons acceptance 


1% (both v and v) 


0.83% (u), 0.22% (P) 


Int. beam v e /v e contamination 


80% (both v and v) 


26% (y), 18% {y) 


Signal normalization error 


5% 


5% 


Background normalization error 


5% 


10% 



Table 1: Detector properties of LArTPC and TASD. 



2.2 Liquid Argon TPC 

Here we consider a 10 kt LArTPC constructed close to NOz^A. LArTPC has excellent 
particle identification and we assume a signal efficiency of 80% for compared to 45% for 
TASD. The energy resolution and background rejection for the LArTPC and TASD are 
comparable. In our simulations, we have taken the efficiencies and migration matrices of 
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LArTPC from [44]. It is expected, of course, that such a detector will come on line much 
later than NOf A. In considering NOi^A + LArTPC, we assume equal 6 years v and v runs 
for NOi/A and equal 3 years v and v runs for LArTPC. The cross-sections for LArTPC 
are slightly different from those of TASD. To obtain the LArTPC cross-sections, we have 
scaled the inclusive charged current (CC) cross sections by 1.06 (0.94) for the v {v) case 
compared to those for water [46, 47]. 

2.3 T2K 

T2K uses the 50 kt Super-Kamiokande water Cerenkov detector (fiducial volume 22.5 kt) 
as the end detector for the neutrino beam from J-PARC. The detector is at a distance of 
295 km from the source at an off-axis angle of 2.5° [16]. The neutrino flux is again peaked 
sharply at the first oscillation of 0.7 GeV. The experiment is scheduled to run for 5 years 
in the neutrino mode with a power of 0.75 MW. Because of the low energy of the peak flux, 
the NC backgrounds are small and they can be rejected based on energy cut. The signal 
efficiency is 87%. The background information and other details are taken from [32, 48]. 

3 Event Rates and Spectrum 

We use GLoBES [49, 50] software to simulate the data for various experiments. For the 
atmospheric/accelerator neutrino parameters, we take the following central (true) values: 

|Am^ ff | = 2.4 • 10 -3 eV 2 , sin 2 26 23 = 1.0, (3.1) 

where Am 2 ff is the effective mass-squared difference measured by the accelerator experi- 
ments in — > disappearance channel [51]. It is related to the AtoJj (larger) and Am^ 
(smaller) mass-square differences through the expression [52] 

Am 2 ff = Am^j — Am 2 ,! (cos 2 6\2 — cos <5cp sin #13 sin 26\2 tan #23), (3-2) 

where Am^ = m 2 . — to 2 . The best fit for #23 is taken from atmospheric neutrino data [53]. 
For #13, we take the global best fit sin 2 20\3 = 0.101 [12] 1 . The uncertainties in the above 
parameters are taken to be cr(sin 2 6>i 3 ) = 13% [12], a(|ATO 2 ff |) = 4% and cr(sin 2 26*23) = 2% 
[16]. We take the solar parameters to be [12] 

Ato^ = 7.62 • 10~ 5 eV 2 , sin 2 9 l2 = 0.32. (3.3) 

We keep these parameters to be fixed throughout the calculation because varying them will 
have negligible effect. We also take the Earth matter density to be a constant 2.8 gm/cc 
because the variations and the uncertainties in density can be neglected for the baselines 
we consider. 

1 The best fit value for sin 2 2#i3 and its uncertainty are taken from the 2nd version of [12], dated 21 May 
2012. The values quoted in the introduction are taken from the 3rd version of the same paper, dated 13 
Aug 2012, by which time our calculations were finished. The 3rd version includes the data presented at the 
Neutrino 2012 conference at Kyoto and hence the values have changed a little. We have checked that our 
conclusions are not affected by these small changes. 
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14 kt TASD @ NOvA Site (NH) 10 kt LAr @ NOvA Site (NH) 




Reconstructed neutrino energy [GeV] Reconstructed neutrino energy [GeV] 

Figure 1: (colour online) Left panel portrays the expected signal and background event rates including 
the efficiency and background rejection capabilities in the v e appearance channel for sin 2 29n = 0.101 and 
Sep = 0°, as a function of the reconstructed neutrino energy for the off-axis NO^A site with 14 kt TASD 
detector running for 3 years exposed to 0.7 MW NuMI beam. Right panel shows the same for a 10 kt 
LArTPC detector running at off-axis NOz^A site for 3 years exposed to 0.7 MW NuMI beam. A normal 
hierarchy has been assumed. In all the panels, the blue dot-dashed and the orange dotted vertical lines 
display the locations of the first and second oscillation maxima. 



The signal and different background event spectra are shown in figure (1) for NOz/A (left 
panel) and for LArTPC (right panel). Even though the mass of the LArTPC detector is 
smaller, the event numbers in the two cases are comparable due to the higher signal ac- 
ceptance. The NC background spectrum is shifted to lower energies due to the use of 
migration matrices and can easily be suppressed by an energy cut. 

In Table-2, we list the number of signal and background events for the old NOz^A detec- 
tor, new reoptimized NOz^A detector and for a 10 kt LArTPC, for both z/„ — > v e appearance 
channel and — > disappearance channel. In the calculation of the event numbers for the 
disappearance channel, we match the numbers with the plot given in NOvA website [54]. 
The background events in the appearance channel arise from three possible sources: (a) 
intrinsic v e /v e content of the beam, (b) mis-identified muons and (c) NC reactions. In 
this study we have taken into account the v e contamination in the beam, which is a 
significant source of the background. In the disappearance channel backgrounds mainly 
arise from the Vu (z^) contamination in (v^) beam with a small number coming from 
NC reactions. The signal events are shown for both neutrino and anti- neutrino runs and 
for both normal and inverted hierarchies. 

We see from Table-2 that both the signal and background event rates are larger for new 
NOz^A compared to old NOz/A. The rise in the signal is only 50% whereas the background 
rises by a factor 3. But the number of signal events increases by 30 whereas the number of 
background events rises by 20. Thus overall, there is a gain in signal relative to the back- 
ground. The background events are fixed and do not change with a change in hierarchy or 
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with a variation of <5cp, whereas the signal events do. Since we are interested in measuring 
these changes in the signal events with hierarchy and Sep, it is more advantageous to have 
a larger signal event sample, even at the expense of a larger background sample. 



Channels 


Old NOuA (15 kt) 


New NOi/A (14 kt) 


LArTPC (10 kt) 


App. 


Signal 


Background 


Signal 


Background 


Signal 


Background 




cc 


(Int+Mis-ID+NC) 


cc 


(Int+Mis-ID+NC) 


cc 


(Int+Mis-ID+NC) 


P Me (NH) 


62 


6+1+4=11 


92 


8+5+19=32 


123 


18+5+7=30 


P Me (IH) 


36 


6+1+4=11 


54 


8+5+19=32 


72 


19+5+7=31 


Ppe(NH) 


26 


6+<l+6=12 


30 


5+<l+10=15 


28 


17+2+2=21 


Ppe(IH) 


34 


5+<l+6=ll 


38 


5+<l+10=15 


36 


14+2+2=18 


Disapp. 


Signal 


Background 


Signal 


Background 


Signal 


Background 




CCQB 


(NC only) 


CCQE 


(NC+Wrong-Sign muon) 


CCQE 


(NC+Wrong-Sign muon) 


Pmm(NH) 


173 


2 


134 


1+6=7 


403 


7+20=27 


Pmm(IH) 


173 


2 


134 


1+6=7 


402 


7+20=27 


Pmm(NH) 


102 


1 


43 


<1+18=18 


136 


2+54=56 


Pw(IH) 


103 


1 


43 


<1+18=18 


137 


2+54=56 



Table 2: Total number of signal and background events for old NO^A, new NO^A and LArTPC, in 
both appearance (z/ M — > v e ) and disappearance (z/ M — > v^) modes. The backgrounds are subdivided into 
three parts: (a) intrinsic beam v e jv e (Int.), (b) mis-identified muons (Mis-ID) and (c) single 7r° events from 
neutral current interactions (NC). In the disappearance mode, wrong-sign muon indicates the background 
coming from (y^) contamination in v^iy^) beam. 




Figure 2: (colour online) Events vs. <5cp- Left panel is for v £ appearance and right panel is for 
appearance. The error bars shown are statistical. Please note the difference in the scale of y-axis in the 
two panels. 

In figure (2), we plot the signal event numbers (along with statistical error bars) as 
a function of 5cp for new NCVA, both for NH and IH. The left panel shows the event 
numbers for three years of running in the neutrino mode and the right panel depicts the 
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event numbers for three years in the anti-neutrino mode. We see that for the combination 
NH and LHP (—180° < 5qp < 0), the neutrino event numbers are much greater than the IH 
event numbers for any <5cp- Similarly for the combination IH and UHP (0 < 5qp < 180°), 
the event numbers are much lower than NH event numbers for any 5qp- Hence NO^A has 
good hierarchy discrimination for these two combinations. For the other two combinations, 
it is possible to have two solutions with different hierarchies for the same event numbers. 
This feature occurs, not only for the total event numbers but also for the spectrum. Thus, 
for unfavourable <5cp half-planes, hierarchy determination is not possible with NO^A. 

4 Numerical Simulation 

We use the minimization of A% 2 to estimate the hierarchy and CP violation sensitivities. 
For hierarchy sensitivity, we first assume NH to be the true hierarchy and we choose a 
true value of 8cp- We compute the NH event spectrum for these assumptions and the 
above true values of neutrino parameters and label it to be data. Then we compute a 
theoretical event spectrum assuming IH and varying the test values of neutrino parameters 
within their ±2<r ranges and <5cp in the full allowed range (—180°, 180°). In doing this 
marginalization, we impose Gaussian priors on the measured neutrino parameters. We 
compute a A% 2 between the event spectra of the data and the theory and demand that 
its minimum value A X 2 min > 2.71 (3.84) for a 90% (95%) C.L. hint of hierarchy. If this 
condition is satisfied, then the hierarchy can be determined at the appropriate confidence 
level, for the assumption of hierarchy being NH and 5qp being the assumed true value. The 
calculation is repeated for various other assumed true values of 5qp- If Ax^m ^ 2.71 (3.84) 
for all values of true Sep, then we can say that the normal hierarchy can be established for 
all possible values of neutrino parameters. The whole calculation then has to be repeated 
for the case of true hierarchy being IH. If Ax^ n ^ 2.71 (3.84) for all true values of Sep, 
then the hierarchy can be determined for all possible values of neutrino parameters. 

In computing the above spectra, it should be noted that the values of Arn^j for NH and 
for IH are different. They are to be calculated from the expression in eq. (3.2) with Am 2 ff 
positive for NH and negative for IH. Since Am|i is always positive, this leads to different 
magnitudes for Am^ for NH and for IH. This difference must be taken into account while 
calculating the NH and IH spectra. Otherwise, there will be a spurious hierarchy sensitivity 
in the disappearance channel. 

As mentioned in the introduction, for a given experiment, we will get two degenerate 
solutions: one with the correct hierarchy and true <5cp and one with the wrong hierarchy 
and a wrong 5cp- It was shown in [31] that two experiments with different baselines, 
with flux peaking at the first oscillation maximum, can pick the correct hierarchy and 5qp 
because the wrong hierarchy solutions for different experiments occur for different values 
of 5qp- But, the statistics of each of the two experiments have to be large enough for a 
clean separation. We will see in the next section that the presently planned runs of T2K 
and NOi^A are not enough for this separation to occur for the full range of 5qp- Therefore, 
for the existence of CP violation, we pose the following question: for what values of true 
Sep can NOzaA. and T2K establish that the CP phase differs from or 180°, independently 
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of the hierarchy? To answer this question, we compute the data spectrum say for NH 
and for a true value of 5qp, as described above. We compute the theoretical spectrum for 
the following four combinations of hierarchy and CP conservation (NH, Sqp = 0), (NH, 
,5 C p = 180°), (IH, 5 C p = 0) and (IH, 5 C p = 180°). We compute the Ax 2 between the data 
and each of these four combinations and choose its minimum value. If Ax^ in ^ 2.71 (3.84), 
then CP violation is established at 90% (95%) C.L. for true hierarchy being NH and for 
Sqp being the assumed true value. We repeat the calculation for the full range of Sep and 
determine the values of true <5cp for which CP violation can be established. As in the case 
of hierarchy determination, these calculations are again repeated for the case where IH is 
the true hierarchy. 

It must be mentioned that this method only establishes that the CP conservation is 
ruled out. But it is possible that the data may not be able to determine the hierarchy. In 
such a situation, the data must be analyzed under the assumption that either of the two 
hierarchies is correct. If the data is analyzed with the assumption of the right hierarchy, 
we get a region of allowed <5cp surrounding the true Sqp ■ If the data is analyzed under the 
assumption of the wrong hierarchy, then the allowed 5qp region may not include the true 
<5cp- The above procedure only requires that, for each of the hierarchies, the allowed Sep 
region should not include the CP conserving cases <5cp = and 180°. If this condition is 
satisfied, the CP violation in neutrino sector is established independent of the hierarchy. 
We will not be able to know the true value of 5cp without determining the hierarchy [35]. 
But proving the existence of of leptonic CP violation will be an important step forward in 
our understanding of the leptonic masses and mixings. 

5 Results 

In this section, we describe the capabilities of the experimental setups considered in this 
work, for the determination of mass hierarchy and CP violation. 

5.1 Mass Hierarchy Discrimination 

In figure (3), we plot the hierarchy discrimination sensitivity of the old NOzM., the new 
NOz^A and the combined sensitivity of new NOi^A and T2K, as a function of the true value 
of <5cp, in the left (right) panel for NH (IH) as the true hierarchy. We see that the wrong 
hierarchy can be ruled out very effectively for 5qp in the favourable half-plane, which is 
LHP (UHP) for NH (IH). The new event selection criteria of NOz^A make the experiment 
even more effective in ruling out the wrong hierarchy for 5cp in the favourable half-plane. 
In the unfavourable half-plane, both the old and the new criteria are equally ineffective. 
However, the addition of T2K data improves the situation significantly and A% 2 increases 
from to > 2 for all the true values of Sep, thus making it possible to get a 90% C.L. 
hint of hierarchy with some additional data. We have checked that a further increment 
in the exposure of T2K or addition of antineutrino data from T2K does not improve the 
hierarchy sensitivity much. 
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Mass Hierarchy Discovery, NH true Mass Hierarchy Discovery, IH true 




8 CP (true) 5 CP (true) 

Figure 3: (colour online) Mass hierarchy discovery as a function of true value of <5cp- Left (right) panel 
is for NH (IH) as true hierarchy. 



In figure (4), we plot the mass hierarchy discrimination capability of NO^A and a stand 
alone LArTPC with three possible masses: 5 kt, 10 kt and 14 kt. Before the reoptimization 
of NOi^A event selection criteria, it was argued that a 5 kt LArTPC has the same capability 
as of the 15 kt TASD detector, because the signal acceptance of the former was three 
times that of the latter [55]. But with the new event selection criteria the performance 
of NOzaA has dramatically improved. We see from the figure (4) that only a LArTPC of 
mass 10 kt can be as effective as NOi^A. But, once again, it must be noted that none of 
the detectors are effective if the true value of 5qp is in the unfavourable half-plane. 
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Figure 4: (colour online) Mass hierarchy discovery as a function of true value of Sep, for different 
LArTPC detector masses and new NOvA. Left (right) panel is for NH (IH) as true hierarchy. 
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Finally, in figure (5) we plot the combined hierarchy discovery sensitivity of NOf A, 
T2K and a LArTPC (of mass 5 kt and 10 kt) and compare it with the hierarchy sensitivity 
of combined NOi^A and T2K data. We find that the addition of even a 5 kt LArTPC 
leads to a very significant improvement in the sensitivity when values of true <5cp are in 
the unfavourable half-plane. With the addition of a 10 kt LArTPC, we see that close to 
95% C.L. hierarchy discrimination becomes possible for all the values of <5cp- 
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Figure 5: (colour online) Mass hierarchy discovery as a function of true value of Sep- Combined 
performance of new NO^A, T2K and LArTPC. Left (right) panel is for NH (IH) as true hierarchy. 



5.2 CP Violation Discovery 



Reoptimization of the event selection criteria of NOi/A has the most dramatic effect on the 
CP violation discovery potential of the experiment. In figure (6), we plot the sensitivity to 
rule out CP conservation scenarios, as a function of true 5cp in the left (right) panel for 
NH (IH) being the true hierarchy. We notice that, while in the case of old NO^A there is 
no CP violation sensitivity at all at 90% C.L., there is such a sensitivity in new NOz/A, for 
about one third fraction of the favourable half-plane. Addition of T2K data leads to CP 
violation sensitivity for about half the region in both favourable half planes at 90% C.L.. 
It can be shown that, T2K by itself, has no CP violation sensitivity. But, the synergistic 
combination of NOi^A and T2K leads to much better CP violation sensitivity compared to 
the individual capabilities. 
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In figure (7), we present a comparison of new NOz^A with stand alone LArTPC of 
different masses. As in the case of the mass hierarchy discrimination, the performance of 
a 10 kt LArTPC detector is closest to the performance of the new NOz^A. 
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Figure 7: (colour online) CP violation discovery as a function of true value of <5cp, for different LArTPC 
detector masses and new NOfA. Left (right) panel is for NH (IH) as true hierarchy. 



In figure (8), we plot the CP violation discovery potential of the combined data from 
new NOz^A with a (6v + 6i>) run, LArTPC of mass 5 kt and 10 kt with a {2>v + 3z?) run 
and T2K along with that of combined new NOz^A with a (3u + 3z?) run and T2K data. 
The net effect of new NOzM. with (6z/ + 6u) run plus a 10 kt LArTPC with a (3f + 3z>) 
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is equivalent to tripling the data of NOi/A with (3f + 3i^) run. We note that the addition 
of the LArTPC leads to a very significant improvement in the range of true <5cp for which 
CP conservation can be ruled out at 95% C.L.. This range goes up from about half of the 
allowed values to about 60% of the allowed values. 
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Figure 8: (colour online) CP violation discovery as a function of true value of Sep ■ Combined performance 
of new NOi^A, T2K and LArTPC. Left (right) panel is for NH (IH) as true hierarchy. 



A summary of our results is given in Table-3. For different combinations of experiments 
and different exposures, we have shown the fraction of 5qp values for which mass hierarchy 
can be determined/CP violation can be detected. 



Setups 


Fraction of <5cp (true) 


MH 


CPV 


NH true 


IH true 


NH true 


IH true 


NOz^A (3+3) 


0.48 (0.43) 


0.46 (0.41) 


0.16 (0) 


0.21(0.04) 


NOvA (3+3) + T2K (5+0) 


0.55 (0.45) 


0.54 (0.43) 


0.38 (0.11) 


0.49 (0.23) 


NOvA (6+6) + T2K (5+0) + 5 kt LArTPC (3+3) 


1 (0.64) 


1 (0.64) 


0.64 (0.56) 


0.68 (0.61) 


NOvA (6+6) + T2K (5+0) + 10 kt LArTPC (3+3) 


1 (0.71) 


1 (0.73) 


0.67 (0.60) 


0.71 (0.64) 



Table 3: Fractions of true values of Sep for which a discovery is possible for MH and CPV. The numbers 
without (with) parentheses correspond to 90% (95%) C.L. The results are shown for both NH and IH as 
true hierarchy. 



6 Summary and Conclusions 

In light of the recent measurements indicating a moderately large value of sin 2 2#i3, the 
NO/vA collaboration has revisited their background rejection cuts. The relaxed acceptance 
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cuts have resulted in a higher number of signal events. Hence, the variations induced by 
the change in hierarchy or a change in 5cp are magnified and become more easy to detect. 

In this work, we have studied the physics reach of the reoptimized NOz/A (in conjunc- 
tion with T2K) to determine the mass hierarchy and CP violation. With the presently 
planned runs of these two experiments, the above goals can be achieved for less than half 
the allowed values of c>cp- There is considerable enthusiasm in the NOz/A and LBNE col- 
laborations for an additional small liquid argon detector detector at the Ash River site. In 
this study we have also considered the possibility of such a module. 

For favourable values of Sep, NOz^A in its original configuration can determine the 
mass hierarchy by itself. The reoptimized new NOz^A allows us to determine the hierarchy 
with greater confidence in the favourable half-plane of Sep- However, in the unfavourable 
half-plane, the sensitivity is negligible. Addition of data from T2K helps to break the 
parameter degeneracy and increases the Ax 2 significantly in the unfavourable half-plane. 
But the currently planned runs of NOz^A and T2K are not sufficient to raise the A% 2 above 
2.71 and give us a 90% hint of hierarchy for 5cp in the unfavourable half-plane. 

A 10 kt LArTPC is found to be equivalent to new NOz/A in its ability to exclude the 
wrong hierarchy. Thus, addition of such a module to the existing NOz/A detector will help 
increase the statistical significance of the experimental data. We find that for increased 
exposure for NOz^A, the combination of data from NOz/A, LArTPC and T2K (nominal) 
can determine hierarchy at almost 2a level for all the values of <5cp- We have checked that 
increasing the exposure of T2K or adding antineutrino data from T2K does not improve 
the results much. 

Discovering CP violation is more difficult than determining the hierarchy. NOz^A by 
itself can discover CP violation only for a small fraction of the favourable half-plane, and 
only at 90% C.L. Addition of data from T2K causes a remarkable increase, by a factor 
of 2.4, in the range of 5qp for which CP violation can be established. This includes a 
significant part of the unfavourable half-plane also. Once again, we find that a 10 kt 
LArTPC has capabilities similar to NOz^A. The combination of additional data from the 
LArTPC detector with boosted NOz/A and nominal T2K greatly improves the ability. For 
around 60% of the entire <5cp range, CP violation can be discovered at 95% confidence 
level. 

In conclusion, we find that after reoptimizing NOz^A for large #13, the sensitivity of the 
experiment to mass hierarchy and CP violation is increased. Adding data from T2K breaks 
the hierarchy-<5cp degeneracy but is not enough to determine hierarchy for unfavourable Sqp 
values. Additional data from a LArTPC will lead to a significant boost in both hierarchy 
and CP violation sensitivities when combined with NOz/A and T2K. 
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